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1^ . Abstract 

' We discuss theoretical uncertainties of the distribution in the cosine of the W po- 

lar angle projected into a measurement of the anomalous triple gauge-boson coupling 
^ . A = = at LEP2 energies for the tandem of the Monte Carlo event genera- 

^ \ tors KoralW and YFSWW3 and for the Monte Carlo event generator RacoonWW. 

Exploiting numerical results of these programs and cross-checks with experimental 
fitting procedures, we estimate that the theoretical uncertainty of the value of A due 
to electroweak corrections, as obtained at LEP2 with the help of these programs, 
is ~ 0.005, about half of the expected experimental error for the combined LEP2 
experiments (~ 0.010). We use certain idealized event selections; however, we argue 
that these results are valid for realistic LEP2 measurements. 

To he submitted to Physics Letters B 
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A distribution of the cosine of the W boson production angle 6w is one of the main 
observables in the VF-pair production process measured by the LEP2 experiments. It 
is sensitive to the triple gauge-boson couplings (TGCs) WWV, with = Z, 7 
particularly to the C- and P-conserving anomalous couplings A^, 0. In this work we 
present our estimate of the theoretical uncertainties (TU) related to electroweak (EW) 
corrections in the measurement of the anomalous TGC A = A-y = A^ in the VT-pair 
production at LEP2 energies^. It is important to note that we treat the anomalous TGC 
as a small new-physics effect beyond the consistent Standard-Model prediction, i.e. we 
introduce the anomalous TGC in the lowest-order matrix element, which is dressed by 
initial-state radiation, while genuine weak corrections are unaffected by the anomalous 
TGC. We concentrate on A in order to keep our analysis as simple as possible and because 
the value of A fitted to experimental data depends more strongly on the shape of the 9w 
distribution than on its total normalization. In fact, while A is mainly sensitive to the 
shape of the W polar angle distribution, it is not particularly sensitive to other single 
distributions. Using only this variable, we lose only 30% in sensitivity On the other 
hand, the analysis of other TGCs would require using more observables. 

Our results were obtained using the Monte Carlo (MC) event generators YFSWW3 ||^ 
g] and KoralW as well as RacoonWW In YFSWW3 and KoralW, 



the anomalous TGCs are implemented according to the most general parametrization 
of ref. and also according to two simplified parametrizations given in ref. 0. In 
Racoon WW0, the anomalous TGCs involving W bosons are implemented according to 
the most general parametrization of ref. |2T| in the form given in ref. 0, while those 
involving only neutral gauge bosons are implemented according to refs. p2|,^ . 



In the following, we consider the semi-leptonic process e+e~ — > M(i/i~z/^, which belongs 
to the so-called CCll class. The corresponding Feynman diagrams constitute a gauge- 
invariant subset of all 4-fermion final-state processes (see e.g. ref. [^] for more details). 
We shall study only the leptonic W polar angle, i.e. the one reconstructed from the four- 
momenta of the fi~ and z/^ (in the actual experiments, the neutrino four-momentum is 
reconstructed from the constrained kinematical fit, see e.g. refs. [0-0]). The input param- 



eters are the same as in the LEP2 MC Workshop studies performed in 2000 |25|. All the 
results in this work are given for the centre-of-mass energy Ecm = 200 GeV. The results 
from YFSWW3 are for the leading-pole approximation (LPA) for the double-resonant 
WW production and decay of the type advocated in ref. [^], which is implemented in 
YFSWW3 as an option LPA^ (see ref. [|12| for more details). The results from RacoonWW 
are based on the double-pole approximation (DPA) as worked out in ref. |TB[. In all our 
numerical exercises we use A = A^ = Xz, and all other anomalous couplings are set to 
zero. 

In our analysis we use two different fitting procedures. In the first part of the paper 
we do direct fits to the cos 6w distributions with the help of the semi-analytical program 
KorWan |p!3|,p!4|,|27[]. Since in KorWan one cannot apply any experimental-like cuts, it can 



-'^A similar analysis for the W mass measurement was done in ref. 0. 

^The anomalous TGCs are implemented in RacoonWW 1.2, which was released recently (see 
ittp : //Itpth . web . psi . ch/racoonww/racoonww . html) . 
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be used, in principle, only for fitting some idealistic distributions (without cuts), for which 
good fits (with low x^) can be expected. We, however, decided to apply the KorWan fits 
also to more experimental-like distributions. The results of these fits should be considered 
only as an independent cross-check for the main fitting procedure, the so-called Monte 
Carlo parametric fit, presented in the second part of the paper. In this procedure, contrary 
to the KorWan fits, any event selection criteria can be applied. As such, this fitting 
procedure can be used for experimental data, which is not the case for the KorWan- 
based method. Note that the results of the KorWan fits depend on the Monte Carlo 
errors, and thus on the models used for the Monte Carlo integration. Since the models 
implemented in KorWan and YFSWW3 are similar, fits of the YFSWW3 distributions 
yield reasonable results for the shifts in A. However, since the models implemented in 
KorWan and Racoon WW are too different to make reasonable fits possible, we do not fit 
the Racoon WW distributions in the first part of the paper. 




Figure 1: Introductory exercise with YFSWW3, see description in text. 

In the first preparatory step, we construct a simple fitting procedure for cos 6*^, where 
the superscript * means that 9w is evaluated in the WW rest frameQ. We "calibrate" the 
fitting procedure for the cos 6*^ distribution by using it with the MC data in which we 
switch on/off the same effects as in the fitting function (FF), typically the initial-state 
radiation (ISR) and the non-factorizable (NF) corrections, in order to see whether we get 
agreement in the case of the same effect in the MC data and in the FF. The FF is taken, in 
all cases, from the semi-analytical program KorWan [0,|14|,^0. In the YFSWW3 MC and 
in the fitting function the NF corrections are implemented in the inclusive approximation 

''More precisely, 6^ is defined as the angle between the directions of the in the WW rest frame 
and the e~ beam in the laboratory frame. 

^The relevant distribution will be available in the next release of KorWan/KoralW. 
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(denoted by INF in the following) of the so-called screened Coulomb ansatz by Chapovsky 
and Khoze which is an approximation of the full calculation of the non-factorizable 



corrections [2S-32]. Here, for completeness, we note that in RacoonWW non-factorizable 
corrections are included beyond the inclusive approximation (see Ref. In partic- 

ular, the real-photonic non-factorizable corrections (as well as final-state radiation) are 
contained in the full e~^e~ — 4/ + 7 matrix elements which are used. 

One immediate profit of this introductory exercise is that we get an estimate of the 
size of the ISR and NF effects on the fitted (measured) A in Fig. |l], which is determined 
from a two-parameter fit to the cos 6^ distribution. In addition to A we also fit the overall 
normalization, which is necessary in the presence of cuts in the MC data and/or when 
the MC model does not correspond exactly to the KorWan model. As a consequence, the 
fitted values of A depend on the normalization and on the shape of the distribution. The 
results of the first exercise, for the BARE47r event selection (no photon recombination) 
and without any cuts (the subscript Atc means the full solid-angle coverage), are shown in 
Fig. |I[ Let us explain briefly the notation: Born denotes the CC03 Born-level results, ISR 
the ones including the 0{a^) LL YFS exponentiation for the ISR as well as the standard 
Coulomb correction [^, INF the above plus the INF correction, and Best denotes the 
best predictions from YFSWW3, i.e. all of the above plus the 0{a^) EW non-leading 
(NL) correction^. In all the MC simulations the input value of A was 0. The errors of 
the fitted values due to the statistical errors of the MC results are in all cases < 0.001. 

Let us summarize the observations resulting from Fig. |^: 

• The fitted values of A agree with the input ones, i.e. A = 0, in the cases when the 
FF corresponds to the model used in the MC simulations. 

• The effects of the ISR are huge, ~ 0.07. 

• The effects of the INF corrections are very small, ~ 0.0015. This is to be expected 
since the non-factorizable corrections are strongly suppressed once the invariant 
masses of the W bosons are integrated over. 

• The effects of the NL corrections are sizeable, ~ 0.008. 

Before we go to the next exercise, let us briefly describe the acceptances/cuts that 
were used in the MC simulations for the following calculations: 

1. All photons within a cone of 5° around the beams were treated as invisible, i.e. their 
momenta were discarded when calculating angles, energies, and invariant masses. 

2. The invariant mass of a visible photon with each charged final-state fermion, Mj^^, 
was calculated, and the minimum value MT^*" was found. If M^^*" < Mr^r or if 

' Jch Jch ^^"^ 

the photon energy < 1 GeV, the photon was combined with the corresponding 
fermion, i.e. the photon four-momentum was added to the fermion four-momentum 



^The 0{a^) EW corrections for the WW production stage in YFSWW3 are based on refs. [p4lp5[ 
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Fits of A, to cos distribution for CAL05 
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Figure 2: Results for calorimetric-type acceptances from YFSWW3. 



and the photon was discarded. In YFSWW3 and KoralW, this was repeated for all 
visible photons, while in Racoon WW there is only one visible photon. 

In our numerical tests we used two values of the recombination cut: 



5GeV: 
25 GeV: 



CAL05, 
CAL025. 



Let us remark that we have changed here the labelling of these recombination cuts 
from the slightly misleading bare and calo names used in Ref. ||2^. They correspond 
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to our CAL05 and CAL025, respectively. This change allows us to reserve the 
BARE name for a "truly bare final fermion" event selection (without any recombi- 
nation) . 

3. Finally, we required that the polar angle of any charged final-state fermion with 
respect to the beams be 6'/^^ > 10°. 

In the next exercise, presented in Fig. 0, we examine similar effects as in Fig. |1], but 
now for the calorimetric-type acceptances/cuts: CAL05 and CAL025. 
The following observations resulting from Fig. § can be made: 

• For the calorimetric event selections, the fitted values of A are shifted by ~ —0.02 
with respect to the input ones in the cases when the FF corresponds to the model 
used in the MC simulations. This can be explained by the fact that the FF of 
KorWan is for the full angular acceptance while the MC results were obtained with 
a cut on the charged final-state fermions. 

• The size of the ISR is the same as for the BARE event selection, cf. Fig. ^ 

• The shift of A due to the NL corrections increased from ~ 0.008 for BARE to ~ 0.011 
for both CAL05 and CAL025. 

• The fitted values of A for CAL05 are slightly different (by ~ 0.002) from the corre- 
sponding ones for CAL025, but the differences AA between various models are the 
same. 

From the fits to KorWan, we can conclude that the effects of the NL corrections on A 
are sizeable, of the order of the expected experimental precision of the final LEP2 data 
analysis. Thus, they need to be accounted for in the experimental measurements of the 
TGCs. 

In the second, main part of our study, we follow an alternative fitting strategy for A in 
which we do not exploit the angular distribution (fitting function) of KorWan, but instead 
we rely entirely on the MC results, more precisely, on a simple polynomial parametrization 
of the normalized angular distribution D{cos9w, ^) = ^ dctsew ^""^^ ^vk, A) as a function 
of cosOw and A, determined from running the MC for several values of A. Since the 
distribution D is normalized, the fitted values of A depend only on the shape. Let us call 
this method a "Monte Carlo parametric fit" , or the MPF in short, and the MPF fitting 
function - the MPFF. Once such a MPF with the MPFF is established, then in principle 
it could be used to quantify a deviation of A by fitting the MPFF to the experimental 
angular distribution. Our aim is rather to quantify various components of TU of the SM 
predictions of the MC tandem KoralW and YFSWW3 and the MC RacoonWW directly in 
terms of A (similarly as in the previous method based on KorWan). The main advantage 
of the MPF is that this can be applied for an arbitrary event selection and any definition 
of 9w, while any semi-analytical program like KorWan has a strongly restricted choice 
of the angle definition and kinematic cuts. The MPF method is feasible and not very 
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difficult in practice because the distribution D{cos6w, A) is a very smooth function of its 
two arguments for LEP2 energies^. Another advantage of the MPF will unfold in the 
following - in fact, we shall be able to get closer to the fitting procedure of the real LEP2 
experiments. 

Here, we include also the results from RacoonWW. The comparison of YFSWW3 and 
Racoon WW is very interesting, because the two calculations differ almost in every aspect 
of the implementation of the ISR, final-state radiation (FSR), NL, and NF corrections. 
The "ISR" prediction of RacoonWW, which also contains the off-shell Coulomb singular- 



ity ||32|,^,^, is exclusively based on the collinear structure-function approach, i.e. all 
generated photons are collinear to the beams and are thus treated as invisible. In contrast, 
YFSWW3 generates photons with finite transverse momenta also in the "ISR" version. 
For the "Best" predictions of RacoonWW, photon radiation in 0{a) is included via the 
full matrix elements, and the main differences to YFSWW3 come from two sources. First, 
RacoonWW does not contain ISR corrections to the non-collinear one-photon emission 
contribution, in contrast to YFSWW3. Then, while the FSR is generated by PHO- 



TOS [§8| in YFSWW3, it is included via the explicit e'^e~ ^ 4/ + 7 matrix elements in 
RacoonWW. 

How is the MPF realized in practice? We use the following 9-parameter MPFF formula 

D(cosO 22 2 

p(cos ew,X)= ( = V aiX' + cos Ow V b.X' + cos' Ow (1) 

^ w y 1=0 j=0 

where parameters Oj, hi and Cj, = 0, 1, 2, are determined by fitting the p(cos6'v(/. A) distri- 
bution obtained from YFSWW3 or RacoonWW for the three values A = —0.2, 0.0, 0.2 []. 
We have checked that the fitted MPFF: (a) reproduces the MC result of p(cos6'ty. A) for 
any A G (-0.2, 0.2) within the fit error of ~ 0.001, and (b) when the MPFF is used to fit 
p(cos6'iy,A) generated by the MC for any A G (—0.2,0.2) the fitted value agrees, within 
the fit error, with the input A used in the MC. Moreover, within fit errors we find ao = 1 
and 60 = Co = as required. Obviously, the other coefficients Oj, hi and q are different for 
every kind of the event selection and angle definition, and the MPF is always a two-step 
procedure: ffist we determine the MPFF using the MC, and then we apply it to the results 
of another MC run (we do not fit experimental data). All the results in the following are 
obtained for the cos 6*^^ distributions, i.e. from now on we use = 

We are now ready to quantify various effects in the cos Q\^^ distribution in terms of A 
using the MPF procedure. To this end, we define pi;'£)(cos6'vi/) = i5i(cos^^, 0)/D2(cos6'ty, 0) 
where D\ and D2 are the normalized angular distributions (for A = 0), whose difference 
we want to quantify and ffi this distribution to the MPFF p(cos^vi/. A) to obtain AA. In 
Table |l|, columns 3-5, we show AA due to the 0[ol) EW NL corrections: in YFSWW3 
denoted by "Y: Best-ISR", in RacoonWW denoted by "R: Best-ISR", and due to the 

^The MPF procedure would be less practical at high energies or even at LEP2 for fitting Mw using 
an effective W mass distribution. 

^The actual fit of a^, hi and a is done in two steps, first for each fixed A — —0.2, 0.0, 0.2, and next 
for the A dependence. We have checked that the results do not change if the fit is done in one step. 
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Fitting procedure 


Fitted data 


Fitting function 


Accept. 


Y: Best-ISR 


R: Best-ISR 


Best: R-Y 


Accept. 


1. KorWan 

2. KorWan 


BARE 
BARE 


0.0114 (6) 
0.0115(6) 


— 


— 


CAL05 
CAL025 


3. MPF: Y-ISR 

4. MPF: Y-ISR 


CAL05 
CAL05 


0.0112(7) 
0.0115(7) 


0.0097(8) 
0.0161 (8) 


0.0008 (9) 
0.0008 (9) 


CAL05 
CAL025 


5. MPF: R-ISR 

6. MPF: R-ISR 


LALUo 
CAL05 


0.0112 (7j 
0.0115(7) 


0.0097 (8j 
0.0161(8) 


0.0008 (9j 
0.0008(10) 


bAL05 
CAL025 


7. MPF: Y-Best 

8. MPF: Y-Best 


CAL05 
CAL05 


0.0113(7) 
0.0116(7) 


0.0098 (8) 
0.0162 (8) 


0.0008(10) 
0.0008 (9) 


CAL05 
CAL025 


9. MPF: R-Best 

10. MPF: R-Best 


CAL05 
CAL05 


0.0110(7) 
0.0113(7) 


0.0096 (8) 
0.0158 (8) 


0.0007(9) 
0.0008 (9) 


CAL05 
CAL025 


11. MPF: KoralW 

12. MPF: KoralW 


ALEPH 
ALEPH 


0.0118(7) 
0.0122(7) 


0.0103 (9) 
0.0172 (9) 


0.0008(10) 
0.0009(10) 


CAL05 
CAL025 



Table 1: The shifts on A from various fits. We use the abbreviations: Y=YFSWW3 and 
R=RacoonWW. The numbers in parentheses are the fit errors corresponding to the last digits of 
the results. 

difference Racoon WW-YFSWW3 in their "Best" modes, denoted by "Best: R-Y". In 
the first two rows of Table 0, we include for the purpose of "backward compatibility" 
AA obtained by using the fitting procedure employing the FF of KorWan (constructed 
originally for BARE4^). Rows 3-12 in Table ^show the results of the MPF procedure. As 
already stressed, the type of the MPF is defined by the variant of the MC and the type 
of the event selection, which are indicated in the first and second columns. In rows 3 and 
4 of Table ^ we use the MPFF determined from the MC run of "Y-ISR", i.e. the MPFF 
was constructed using the results from YFSWW3 in the "ISR" mode (described above). 
The similar MPFF from RacoonWW, denoted by "R-ISR" was used to obtain the results 
in rows 5 and 6. In rows 7 and 8 "Y-Best" denotes the MPFF constructed from the run of 
YFSWW3 in the "Best" mode, and the similar one for RacoonWW corresponds to rows 9 
and 10 ("R-Best"). The type of the event selection (CAL05 or CAL025) in the MC data 
is indicated in the last column. We use the MPFFs with the coefficients for CAL05 as it 
gives practically the same results for the fitted A as the ones for CAL025. We see that 
all the results of the above five MPFs are consistent with each other, and for YFSWW3 
they also agree quite well with the results of the fits employing the FF of KorWan. 

The results in Table |I| show that the NL corrections influence A considerably; they shift 
A by —0.0096 to —0.0172 with respect to the ISR-level predictions, which is comparable 
to the expected final experimental precision on this TGC at LEP2. The large differences 
in the shifts of A for "R: Best-ISR" between CAL05 and CAL025 originate from the fact 
that the "ISR" prediction of RacoonWW includes only coUinear photons and thus is not 
influenced by photon recombination, in contrast to the "Best" prediction of RacoonWW 
and the "ISR" and "Best" predictions of YFSWW. The results of Table demonstrate 
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that the NL corrections have to be included in the respective data analyses. On the other 
hand, the differences between the "Best" predictions of Racoon WW and YFSWW3 are 
very small, < 0.001. This is important as the two programs differ in many aspects of 
the implementation of various effects. The good agreement between these programs is an 
indication that the effects on A due to the used approximations and the missing higher- 
order corrections should be small. This will be further investigated in the following. 

Last but not least, the MPFF can also be defined for the true LEP2 acceptance. In 
rows 11 and 12 we exploit MPFFs that use Oj, bi and q determined for the real-life ALEPH 
acceptance, which are denoted as "ALEPH". The ALEPH selection and reconstruction 
efficiencies for jets and leptons are described in refs. [|l],^ (typically, the lepton acceptance 
reaches |cos^| = 0.95). Again the values of the fitted AA are about the same as in the 
previous exercises for the "academic" event selections. This important cross-check makes 
us confident that our estimates of the TU of A are indeed relevant to the real-life LEP2 
measurements. 



Effect 


Acceptance 


AA 


1. Best - ISR 


BARE4,r 

CAL054^ 


0.0108(7) 
0.0110(7) 


2. ISR3 - ISR2 


BARE4,r 

CAL054^ 


0.0001 (2) 
0.0001 (2) 


3. FSR2 — FSRi 


BARE4,r 

CAL054^ 


0.0001 (3) 
0.0001 (3) 


4. 4/-background corr. (Born) 


CAL05 
CAL025 


0.0021 (3) 
0.0021 (3) 


5. 4/-background corr. (with ISR) 


CAL05 
CAL025 


0.0005 (3) 
0.0005 (3) 


6. EW-corr. scheme: {B) — (A) 


CAL05 
CAL025 


0.0006 (9) 
0.0006 (9) 


7. LFAb - LFAa 


CAL05 
CAL025 


0.0017(9) 
0.0018(9) 



Table 2: The shifts on A from various effects, obtained with YFSWW3 and KoralW. See the text 
for more details. 

In Table 0, we present the results of the A-shifts due to switching on/off various 
effects in YFSWW3 and KoralW. They were obtained with the MPFF constructed with 
YFSWW3-Best for CAL05. The subscript in in the acceptance denotes the full sohd- 
angle coverage. In row 1, we show the results of the MPF for the NL correction in 
YFSWW3, but for acceptances slightly different from the ones in Table ^ namely for 
BARE47r (no photon recombination) and CAL0547r (similar to CAL05, but without cuts 
1 and 3). These results agree very well with the corresponding results in Table |l] for 
CAL05 and CAL025 (cf. column 3). This shows that the value of A is not very sensitive 
to the choice of cuts and acceptances. Rows 2 and 3 contain the results of switching off 
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the 3rd order ISR and the 2nd order FSR corrections, respectively. Both these effects 
are neghgible at LEP2 in terms of the A-shifts. The 4/-background corrections^ are 
investigated in rows 4 and 5. While a 4/-background contribution at the Born level 
induces a shift on A of ~ 0.002, it gives a negligible effect when it is combined with the 
ISR. This is because the ISR affects the cos 6'^^ distribution in a way opposite to that 
of the 4/-background, which leads to large cancellations of the latter effect. Then we 
investigate the uncertainties of the NL corrections due to different schemes for the EW 
effective couplings (they account for some parts of higher-order EW corrections): the 



so-called schemes (A) and (B) in YFSWW3 |jT2[ (row 6), and due to the different LPA 



approaches: LPA^ versus LPA;, (row 7). While the effects on A of the differences 
between schemes (A) and (B) are negligible (below the fit errors), the ones due to the 
LPA variation are ~ 0.002 - this can be regarded as the LPA uncertainty in YFSWW3. 



Effect 


Acceptance 


AA 


1. Best - ISR 


CAL05 
CAL025 


0.0096 (8) 
0.0158 (8) 


2. Off-shell Coulomb effect 


CAL05 
CAL025 


0.0001 (10) 
0.0001 (10) 


3. 4/-background corr. (Born) 


CAL05 
CAL025 


0.0029(10) 
0.0029(10) 


4. 4/-background corr. (with ISR) 


CAL05 
CAL025 


0.0008(10) 
0.0008(10) 


5. On-shell projection 


CAL05 
CAL025 


0.0003(10) 
0.0003(10) 


6. DPA definition 


CAL05 
CAL025 


0.0005(10) 
0.0005(10) 



Table 3: The shifts on A from various effects, obtained with RacoonWW. See the text for more 
details. 

In Table ^, we show results of the A-shifts due to switching on/off various effects in 
RacoonWW, as obtained with the MPFF constructed with RacoonWW-Best for CAL05. 
In row 1, we give the shifts between the "Best" and "ISR" modes of RacoonWW. These 
shifts are due to the NL electroweak corrections, including, in contrast to the numbers 
from YFSWW3 in Table ^, also the effects of non-collinear photon emission. Row 2 
exhibits the results of switching off the off-shellness of the Coulomb singularity. This 
effect is negligible, since the energy is not in the threshold region and the Coulomb 
singularity affects mainly the normalization and thus cancels in D and p. The effect of 
the 4/-background diagrams without and with the corresponding ISR (rows 3 and 4) is 
consistent with the corresponding results in Table |^. Finally, we study the uncertainties 

®The complete Born-level 4/ matrix element in KoralW was generated with the help of the GRACE2 
package |4C| ]. 
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of the NL corrections due to different definitions of the DPA in RacoonWW (rows 5 and 
6) [jl8| . These effects are smaller than the corresponding effects seen in row 7 of Table 0. 
This is because, in RacoonWW, only the virtual corrections are treated in DPA, while in 
YFSWW3 also the real corrections feel the LPA. 



ALEPH MPFF 


Fitted data 


Channel 


Acceptance 


Y: Best-ISR 


R: Best-ISR 


Best: R-Y 


Acceptance 




TRUE 


0.0118(7) 
0.0121 (7) 


0.0102 (9) 
0.0170 (9) 


0.0008(10) 
0.0009(10) 


CAL05 
CAL025 


RECO 


0.0118(7) 
0.0122(7) 


0.0103 (9) 
0.0172 (9) 


0.0008(10) 
0.0009(10) 


CAL05 
CAL025 


eveqq 


TRUE 


0.0119(7) 
0.0122(7) 


0.0103 (9) 
0.0172 (9) 


0.0008(10) 
0.0009(10) 


CAL05 
CAL025 


RECO 


0.0119(7) 
0.0123(7) 


0.0103 (9) 
0.0172 (9) 


0.0008(10) 
0.0009(10) 


CAL05 
CAL025 


TUrqq 


TRUE 


0.0115(7) 
0.0118(7) 


0.0100 (8) 
0.0166 (8) 


0.0008(10) 
0.0009(10) 


CAL05 
CAL025 


RECO 


0.0107(6) 
0.0109 (6) 


0.0091 (8) 
0.0152 (8) 


0.0007(9) 
0.0008 (9) 


CAL05reco 
CAL025reco 


qqqq 


TRUE 


0.0118(7) 
0.0120(7) 


0.0102 (9) 
0.0169 (9) 


0.0008(10) 
0.0009(10) 


CAL05 
CAL025 


RECO 


0.0094 (6) 
0.0096 (6) 


0.0081 (7) 
0.0132(7) 


0.0007(8) 
0.0008 (8) 


CAL05reco 
CAL025reco 



Table 4: The results for A shifts for the fits using the ALEPH fitting functions obtained at the parton 
level (TRUE) and with the full detector reconstruction (RECO) for various channels. The subscript 
RECO for some CAL05 and CAL025 acceptances means that the reconstruction effects were also 
included in the corresponding MC data (see the text for details). The fitted data are always for the 
channel fi^i?^ud. We use the abbreviations: Y=YFSWW3 and R=RacoonWW. 

In order to see whether our results can be extended also to other channels, we per- 
formed the following exercise. Using simulated data at the parton level (denoted as 
TRUE) and with the full ALEPH detector reconstruction (denoted as RECO) for the 
channels l^i^^qq, ^^eQQ, t^tQQ and qqqq, we constructed eight MPFFs - two MPFFs per 
channel. Then we fitted these functions to our MC data, the same as in Table ^ for the 
channel ix'v^ud. The results of these fits for all the TRUE-level fitting functions and 
for the ^iv^qq and eveqq RECO-level ones are shown in Table |. As one can see, all the 
TRUE-level results are consistent with each other and agree with the ones for the YF- 
SWW3 and RacoonWW MPFFs (cf. Table |, rows 3-10). The RECO-level results agree 
very well with the TRUE-level ones for the channels fn^^qq and eUeqq. The above results 
show that the experimentally reconstructed distributions of cos 9w for the channels iJiv^qq 
and ez/e^5 are very similar to the parton-level ones and also to each other, which in our 
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exercise resulted in almost identical fitting functions. This leads us to the conclusion that 
our previous findings for the channel fJ^i^^qq can be extended to the channel eUgQQ- 

For the ruT-qq and qqqq channels, the reconstructed distributions differ considerably 
from the TRUE-level ones. We have checked that fitting the RECO-level MPFF from 
the qqqq {Tv^qq) channel to our MC data results in the shifts of A which are by ~ 100% 
(~ 25%) larger than the corresponding ones for the TRUE-level MPFF (these results are 
not shown in Table Thus, the parton-level MC data are not appropriate for estimating 
the shifts of A in the real-life experimental measurements for these channels. In order to 
obtain realistic effects, these data have to be processed through the full detector simula- 
tion. Instead of feeding our MC events into the full ALEPH reconstruction generator, we 
applied so-called transfer matrices to the MC cos9]j^^ distributions. A transfer matrix 
gives the probability for an event generated in a TRUE cos 6y^^ bin to be found in a given 
RECO bin. This takes into account non-diagonal terms induced by the ALEPH detector 
effects (jet resolution, jet pairing, jet charge, and r reconstruction) |jl|. The appropri- 
ate transfer matrices for these channels were constructed using the Monte Carlo samples 



generated within the full ALEPH detector simulation environment ^1[]. Then, we used 
such experimental-like distributions in the MPFs with the RECO-level fitting functions. 
In these fits we used appropriate error (covariance) matrices, taking into account corre- 
lations between the bins and between the transfer-matrix elements. The results for these 
two channels at the RECO-level complete Table ^ (they are denoted by the subscript 
RECO attached to the respective acceptance in the last column). As can be seen, the 
values of the A shifts are close to the corresponding TRUE-level results. We also checked 
that this method did not affect the ^iv^^qq and eveqq channels, as expected. Therefore our 
conclusions for the TU of A in the fi~Ufj_ud channel can be applied also to the rurqq and 
qqqq channels in the realistic experimental analyses. 

From the above numerical exercises and the accompanying discussion, we come to the 
following conclusions: 

• The non-leading EW corrections cause shifts of A = A^ = A^ at the level of 0.01-0.02, 
which is comparable to the combined experimental accuracy for this anomalous TGC 
at LEP2. Thus, these corrections have to be taken into account in the experimental 
analyses. 

• The comparisons between YFSWW3 and RacoonWW and the tests of various 
modes/options in both programs (cf. Tables and ||) allow us to estimate the EW 
theoretical uncertainty in A. We use the largest shifts found in these comparisons 
and apply a safety factor of 2 to account for the ~ 30% sensitivity loss due to the 
single-distribution fit and for possible higher-order effects missing in both programs. 

From this we estimate the EW theoretical uncertainty in A of the MC tandem 
KoralW&YFSWWS and of the MC program RacoonWW to be ~ 0.005 at LEP2 
energies, which is ~ 1/2 of the expected combined experimental error. 
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